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1. Introduction 

The third Fermi Large Area Telescope catalogue (3FGL) of high energy y-ray sources contains 
3033 sources found during 4 years of operation of the instrument [|I|]. Out of these, 1009 objects are 
unidentified as no counterparts at other wavelengths have been associated with them. More sources 
are yet to be revealed since Fermi-LAT has been conducting observations for almost 7 years. 

In this work, we present our finding of a previously undocumented source of a signal in the 
Fermi-LAT data. It emerged during our analysis of a region around the BL Lac lES 0229+200, 
our initial source of interest. Moreover, several other fainter objects were revealed along that. The 
analysis and its results are described in Section The discussion of the features of the found y-ray 
excesses is given in Section 

2. Analysis 

We used the publicly available Pass 7 data from a region of interest (ROI) with a radius of 
15°. The centre of ROI was at the position of lES 0229+200 at RA = 02h 32m 48.6s, Dec = 
+20° 17' 17.5". 

We examined the data using the standard analysis chain of Eermi Science Tools (v9r33p0) 
The chosen time range encompasses events registered between August 4, 2008 and October 31, 
2014. Events with energies between 100 MeV and 100 GeV were selected. A cut on the maximum 
zenith angle of 100° was applied. After the preliminary selections, we proceeded with the binned 
analysis and used the recommended instrument response functions P7REP_SOURCE_V15. 

The contribution from nearby sources to the observed emission in the ROI was accounted 
for by adding their spectral and spatial descriptions to the source model. The diffuse background 
was modelled using the templates gll_iem_v05_revl and iso_source_v0 5_revl for the 
galactic and isotropic diffuse emission, respectively. The model also contained known sources up 
to 25° from the centre of the ROI in order to take into account the point spread function (PSP) 
of Permi-EAT at low energies [^. The list of sources together with their spectral and spatial 
specifications was initially taken from the second Permi-EAT catalogue (2PGE) made after 2 years 
of data taking 0. All spectral parameters were kept fixed for very faint sources and for those 
which lie beyond 7° from the centre of the ROI. Since IPS 0229+200 is not present in the 2PGE 
catalogue, we added another source in the model at the centre of ROI and assumed its spectral 
shape is a simple power law. 

Maximum likelihood estimates (MEE) of relevant parameters were found using the gtlike tool. 
In the preliminary analysis, we kept in the model only those sources, the contributions of which 
to the emission in the ROI were not negligible. In the end, the blazar IPS 0229+200 was firmly 
recognized in the data with the test statistic TS = —2 (log.ifo “ log-5f) = 70. Here, .if and Jfo are 
the likelihood functions for the model with and without the additional source included, respectively. 

2.1 Source identification 

We checked for any additional sources which were not included in our source model but, 
nonetheless, contribute significantly to the observed signal in the investigated region. Por this 

* ittp://fermi.gsfc.nasa.gov/ssc/data/analysis/softwara 
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Figure 1: A map of test statistic values for bins of 0.25° size in a 15° x 15° region centred around the 
position of lES 0229+200. Only sources present in the 2FGL catalogue were assumed to contribute 
to the /-ray signal in the shown region. Vertical bar at the right side represents colour encoding for 
different TS values. Gaussian smoothing was applied to the map. Labelled crosses denote sources 
present in the 3FGL but not in the 2FGL catalogue. White circle marks the unexpected signal 
labelled as S1 unattributable to any catalogue source. 


purpose, we constructed a map of TS values based on the output model file which was obtained 
after the iterations of the gtlike analysis. lES 0229+200 was excluded from this model. We used 
the gttsmap tool which calculates the TS values for individual spatial bins by introducing a putative 
source at a given position. The map covers an area of 15° x 15° and its binning was chosen as 0.25° 
per pixel. 

The TS map is shown in Eig. [I[ The area in the lower left corner of the map is the only spot 
of residual TS values that can be linked to a source from the 2EGL catalogue. We believe that this 
bright area is a relic of 2EGL J0211.2+1050 which is located beyond the boundary of the map. 
Although this source was included in our model, we presume that it appeared in the map due to its 
flaring activity in 2010 and 2011. 

In Eig. H, the presence of lES 0229+200 in the data is visible as the spot in the centre of 
the map. Additional bright areas appear in the TS map which do not correspond to any 2EGL 
counterpart. After the release of the 3EGL catalogue [|^], we compared the positions of these 
hotspots with the coordinates of sources reported in the 3EGL but not in the 2EGL catalogue (green 
crosses in Eig. ||). In the investigated region, most of the emitters from the newer catalogue can 
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Figure 2: A map of test statistic values based on a source model including sources from the 3FGL 
catalogue and the hypothetical object SI. Circles encompass the unknown sources of HE /-rays. 
Green crosses pinpoint the positions of two known BE Eac objects [Q]. Eor further details see 
caption to Eig. [I[ 

be easily paired with the sources which appear in our TS map. Note that the catalogue name of 
lES 0229+200 is 3EGE J0232.8+2016. 

The area labelled as SI in Eig. |T] cannot be attributed to any of the known HE y-ray sources. 
This was confirmed after we rerun the analysis with the new source model made out of all 3EGE 
sources situated within the 25° radius. S1 remained the only clearly distinguishable residual signal. 

In the next step, we included an additional source in our model. Its spectral type was assumed 
to be a power law and the coordinates were read off from the TS map as the approximate position of 
the S1 hotspot. After the MEE of the spectral parameters were found, gtfindsrc tool was employed 
to find the best estimate of the position of the presumed source. The resulting coordinates of S1 
along with the TS value associated with this additional source are given in Table |T| 

In Eig. 1^ a TS map obtained for the source model with SI included is shown. No other 
significant residuals can be recognized around its position. However, four more spots became 
distinguishable at levels of significance sufficient for claims of detection (white circles in Eig. 

We performed the same analysis routine as described above for the four new source candidates. In 
Table |I|, the coordinates and TS values of these objects are listed. 

Einally, the gtlike tool was retried with a model containing all 3EGE sources as well as five ob¬ 
jects S1-S5. All of them were assumed to be point-like with power law shaped spectra. The results 
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Object 

TS 

RA [°] 

Dec [°] 

SI 

90.3 

36.09 

18.53 

S2 

36.3 

32.51 

25.31 

S3 

41.2 

34.13 

23.25 

S4 

31.6 

42.03 

22.54 

S5 

34.0 

40.62 

17.57 


Table 1: Test statistics associated with the detection of five unknown sources and their positions. 


Object 

No 

(ph cm^^ s^^ MeV^*) 

r 

/ (> 100 MeV) 

(ph cm-^ s-*) 

SI 

(4.03 ±0.67) X 10-^° 

2.81 ±0.09 

(2.22 ±0.30) X 10-^ 

S2 

(1.71 ±1.06) X 10-“ 

1.87±0.16 

(1.96 ±0.90) X 10-5 

S3 

(1.40±0.89) X 10-“ 

1.81±0.16 

(1.72 ±0.79) X 10-5 

S4 

(1.62±1.15) X 10-“ 

1.86±0.17 

(1.88±0.99) X 10-5 

S5 

(1.14±0.60) X 10-*° 

2.41±0.18 

(8.09±3.36) X 10-5 

lES 0229+200 (4.80± 1.76) x 10’^^ 

2.04±0.13 

(4.61 ±1.31) X 10-5 


Table 2: Maximum likelihood estimates of spectral parameters of analysed objects. Differential 
energy spectra are given by the power law prescription as F{E) = No {E/EqY^■ Here, No is the 
flux normalization at Eq = 100 MeV and T is the photon index. Integral flux above 100 MeV is 
given in the last column. 


of the analysis for the new objects are given in Table Spectral parameters of lES 0229+200 are 
also stated therein. 


3. Discussion 


Particular attention was paid to the characteristics of the initial subject of investigation, the 
blazar lES 0229+200, and the brightest of the new hotspots, SI. We examine the spectral energy 
distributions of both sources in Section 3.1. Changes of the their fluxes over time are inspected in 
Section 3.2. Discussion of detection of SI is given in Section 3.3. Einally, possible astrophysical 


counterparts to the unknown source candidates are reported in Section 3.4 


3.1 Spectral energy distributions 

Spectral pai'ameters of lES 0229+200 deduced from over 6 years of data (see Table §) are 
in good agreement with those reported by the Eermi-EAT collaboration after 4 years of obser¬ 
vations [||]. Spectral energy distribution of lES 0229+200 is depicted in Eig. The source is 
not detected below 300 MeV. Maximum likelihood estimates of the spectral parameters obtained 
for the data within the range 0.1 — 100 GeV suggest that the source exhibits a hard energy spec¬ 
trum not inconsistent with the power law shape. On closer inspection, we notice that the source 
spectrum is steeper in the range 0.3 — 2 GeV than at higher energies. The high energy part of the 
spectrum is compatible with the results of Vovk et al. (2012) who studied the source in the range 
1-300 GeV0. 
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Figure 3: Spectral energy distributions of lES 0229+200 and SI. Upper limits on the flux at a 95% 
confidence level denoted by arrows are stated in those bins where TS < 1. Solid lines represent 
the spectral shapes given by maximum likelihood estimates of the relevant pai'ameters obtained 
for the data from the whole studied energy range. Grey areas are defined by uncertainties on the 
differential fluxes. The error bars denote the statistical uncertainties at a 1 a confidence level. 


The self-synchrotron Compton (SSC) scenario is commonly employed to explain the origin 
of /-radiation in blazars. The one-zone SSC model [|5|] was used by the VERITAS collaboration 
to describe the available multiwavelength data [^, including the results of previous analysis of 
Eermi-EAT data 0. The GeV-TeV data were well described by a single peak attributed to the 
inversely Compton scattered population of photons (see Eig. 6 in Ref. [^). However, our results 
suggest that a break may be present in the spectrum around the energy of 2 GeV (see Eig. This 
may be a sign of a secondary component in the source /-ray emission. 

S1 is most active at energies below roughly 3 GeV, see Eig. |3^. The spectral analysis revealed 
that the test statistic corresponding to the detection of S1 above 3 GeV is TS < 1. The signal is most 
significant in the energy range where the containment angle of photon directions reconstructed by 
Eermi-EAT is the largest. Thus, it is expectable that the extent of S1 is the largest among ah new 
hotspots. In Eig. ^ the lack of any residual signal in the surroundings of SI (compare with Eig. 
demonstrates that the addition of a hypothetical point-like object at the position of S1 can indeed 
fully account for the unknown signal. 

3.2 Temporal evolution 

Temporal evolution of the integral fluxes of lES 0229+200 and SI is depicted in Eig. ^ The 
whole investigated time range was divided into six equally long periods. Eor each of them, integral 
fluxes were calculated assuming constant spectral indices E as given in Table 

The fit of lES 0229+200 data to a constant yields = 4.4 for 5 degrees of freedom. The 
probability that such or greater value of x'^ statistic is observed under the null hypothesis of stable 
flux is P{x^,5) Ri 0.5. No significant deviations from the constant flux on a time scale of years can 
be recognized in the case of lES 0229+200. 

We note, however, that there has been some evidence for annual variability of this source in 
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Figure 4: Lightcurves of lES 0229+200 and SI. The dashed line represents the fit of the data 
from the direction of lES 0229+200 to a constant. The dash-dotted line depicts the linear fit of SI 
data. Statistical uncertainties at a la confidence level are sfated. 


fhe very high energy (VHE) regime above 300 GeV [^, In particular, the VHE intensity of the 
source was found to be in a higher state in the 2009-2010 season (MJD 55118-55212) [^. 

The y-ray intensity of S1 was changing over time, see Eig. The x^ fit of data from the di¬ 
rection of SI to a constant (dashed line) yields a of 19.3 for 5 degrees of freedom corresponding 
to the y>-value of 1.7 x 10^^. As the flux seems to follow an increasing trend over time, we fitted a 
linear function to the data points, see the dash-dotted line in Eig. ^ Eor the linear fit, we obtained 
a x'^ value of 6.1 for 4 degrees of freedom corresponding to the p-value of 0.2. 

3.3 Detection of SI 

It is intriguing that though SI comes out to be brighter than lES 0229+200 in our analysis, SI 
was not added to the 3EGE catalogue ||^] like the latter source. We inspected this by analysing 4 
years of data, i.e. the time range examined by the Eermi-EAT collaboration for the creation of the 
catalogue. 

We used a source model containing all emitters from the 3EGE catalogue |[T]] to test the alter¬ 
native hypotheses that new sources lES 0229+200 and SI contribute to the signal in the studied 
region. We obtained test statistics TS = 48 and 20 for lES 0229+200 and SI, respectively. Thus, 
when taking into account contributions from all new source candidates in ROI, the test statistic for 
S1 is less than the minimal value of 25 that was set by the Eermi-EAT collaboration for keeping 
the source candidates in the 3EGE catalogue |[I]]. It is also worth noting that the activity of SI 
increased in the following years and the object became easier to recognize, see Eig. and Table ||. 

3.4 Counterparts 

We investigated the possible physical origin of y-ray sources revealed in our analysis. Eor 
that purpose, we looked for astrophysical objects which may be consistent with the positions of the 
suspected sources. We found that two of the hotspots, S3 and S5, are very well compatible with 
the positions of BE Eac objects RXS J02165+2314 and TEX 0239+175 0, respectively. These 
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active galaxies are marked by green crosses in Fig. Q To our knowledge, neither of them has been 
reported to be active in the high energy regime so far. 

4. Conclusions 

We analysed 6.2 years of high energy y-ray data gathered by the Fermi-LAT experiment 
from a region surrounding the BL Lac lES 0229+200. The features of the source were found to 
be consistent with older findings ||I|, We noticed a possible break in the spectrum around the 
energy of 2 GeV. No increase of the source activity in the periods of the VHE enhancement was 
observed. 

We found an unexpected significant signal near lES 0229+200. Additional four fainter un¬ 
known sources were recognized in the studied region. None of the revealed objects have been 
previously reported. Spectral shapes and fluxes of new sources of a y-ray excess were recorded. 
Possible correlation of fwo signal areas wifh BL Lacs RXS J02165+2314 and TEX 0239+175 was 
found. The brighfesf of fhe unknown signals was sfudied in defail. No support for fhe extended 
nafure of fhis candidafe source has been esfablished. Ifs emission is variable and can be described 
by a steep power law. 
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